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ABSTRACT 

Standard f)o = 1 cold dark matter (CDM) needs 0.27 < ag < 0.63 (2o-) to fit the 
observed number of large separation lenses, and the constraint is nearly independent 
of Hq = 100/i~^ km s~^ Mpc~^. This range is strongly inconsistent with the COBE 
estimate of ag = (2.8 ± 0.2)h. Tilting the primordial spectrum oc k^ from ra = 1 
to 0.3 ^ ra ^ 0.7, using an effective Hubble constant of 0.15 ^ T = h ^ 0.30, or 
reducing the matter density to 0.15 ^ ^oh ^ 0.3 either with no cosmological constant 
(Ao = 0) or in a fiat universe with a cosmological constant (fio + Aq = 1) can bring 
the lensing estimate of ag into agreement with the COBE estimates. The models 
and values for ag consistent with both lensing and COBE match the estimates from 
the local number density of clusters and correlation functions. The conclusions are 
insensitive to systematic errors except for the assumption that cluster core radii are 
singular. If clusters with p oc (r^ + s^)~^ have core radii exceeding s = 15h~^a^ kpc 
for a cluster with velocity dispersion a = 10"^ (T3 km s~^ then the estimates are invalid. 
There is, however, a fine tuning problem in making the cluster core radii large enough 
to invalidate the estimates of ag while producing several lenses that do not have central 
or "odd images." The estimated completeness of the current samples of lenses larger 
than 5'.'0 is 20%, because neither quasar surveys nor lens surveys are optimized to find 
this class of lenses. 



Subject headings: gravitational lensing - cosmology: observations - cosmology: theory 
- dark matter - large-scale structure of universe 

1 Introduction 

Narayan & White (1988) pointed out that the standard cold dark matter (CDM) model 
predicts many more gravitational lenses with separations above 5'.'0 than we find. Their 
predictions were based on the Press-Schechter (1974, PS hereafter) model, but more recent 
numerical models by Cen et al. (1994) and Wambsganss et al. (1994) confirm the fun- 
damental result. A singular isothermal sphere with velocity dispersion a = 10"^ (T3 km s~^ 
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produces lenses with an average image separation of A^ = 28'.'8(t|, so the large separation 
lenses explore the number and evolution of clusters and groups. They are, however, a qual- 
itatively different test of cosmogonic models than the local density of clusters (Peebles et 
al. f989, Frenk et al. f990, Bahcall & Cen f992, f993) or correlation functions (Maddox et 
al. f990, Picard f99f, Vogeley et al. f992, Loveday et al. f992). Any massive, collapsed, 
virialized halo will produce lenses, so the test is independent of the luminosity of the lenses 
and unaffected by problems with detecting and counting complete cluster samples locally. 
The probability of lensing peaks at intermediate redshifts, and it goes to zero at low redshift, 
so it is a test of the number density of groups and clusters at z ^ 0.3-0.5 rather than at 
z = 0. This allows lensing to distinguish between scenarios that produce the same number 
of clusters today using different formation histories. 

If we can understand the selection function for large separation lenses in heterogeneous 
quasar catalogs, compute the magnification bias of the sample, and decide which of the 
large separation quasar pairs to call lenses, then we have an important new cosmological 
probe. Only qualitative comparisons can be made without including magnification bias and 
selection effects, because they can change the number of lenses found in any observational 
sample by an order of magnitude. Unfortunately, Narayan & White (1988), Cen et al. 
(1994), and Wambsganss et al. (1994) only made qualitative comparisons between the models 
and the observations. Narayan & White (1988) and Wambsganss et al. (1994) did not 
include selection effects and magnification bias, and Cen et al. (1994) used a crude model 
based on lens surveys. Where Cen et al. (1994) and Wambsganss et al. (1994) pursued 
numerical calculations of cross sections, we will focus on selection effects, magnification bias, 
and quantitative estimates of the number of lenses in various cosmological scenarios. The 
disadvantage of our approach is that, like Narayan & White (1988), we rely on the PS model 
to estimate the number and distribution of lenses. We can, however, see if the PS formalism 
is an accurate method for estimating lens probabilities by comparing to the Wambsganss et 
al. (1994) numerical results when possible and by examining the effects of the systematic 
uncertainties on the conclusions. The advantage of the PS model is that we can rapidly 
survey a large number of cosmological scenarios to examine the sensitivity of the method to 
its parameters and to see which models are constrained by gravitational lensing. 

We know of two confirmed lenses with separations larger than S'.'O (Q0957-|-56I and 
Q20I6+II2), and another four candidates (QII20+0I9=UM 425, QI429-008, QI635+267, 
and Q2345+007). Two additional pairs, PKS 1 145-071 AB (Djorgovski et al. 1987) and 
QI343-|-266AB (Crampton et al. 1988), are rejected as lens candidates even though the 
redshift difference in both pairs is less than Az < O.OOI. PKS 1 145— 071 is rejected because 
one quasar is radio loud and the other is radio quiet (>500:I fiux ratio), and QI343-|-266AB 
is rejected because of gross differences in the spectral lines of the quasars. The properties 
of these eight objects are summarized in Table I. A key distinction between the pairs in 
Table I is whether they were found as part of the original survey that found the quasar, or 
whether they were found in a lens survey examining known quasars to see if they are lensed. 
Four of the eight objects in Table I were found in the original quasar survey (Q0957-|-56I, 
QI343-|-266, QI635-I-267, and Q2345-I-007) and four were found as part of a search for lensed 
images (QII20+0I9, PKS II45-07I, QI429-008, and Q20I6+II2). QII20+0I9 and PKS 



1145—071 were found in a survey by Djorgovski & Meylan (1989), Q1429— 008 in a survey 
by Webster et al. (1988), and Q2016+112 in the MG survey (Burke et al 1992). In the 
first two cases the survey fists and selection functions are unpubfisfied, so we cannot build a 
tfieoretical modefi For Q2016+112 we do not fiave tfie necessary information on tfie redsfiift 
and radio fiux distributions of tfie MG sources to make a tfieoretical model. For comparison, 
all tfie galaxy scale lenses but one (PGl 115+080, Weymann et al. 1980) were found in lens 
surveys. 

Tfie key to drawing quantitative conclusions is §2, wfiere we develop a selection effects 
model for finding lenses in a fieterogeneous quasar like tfie Hewitt-Burbidge (1993, HB93 
fiereafter) catalog. More importantly, we sfiow tfiat it is a valid selection effects model for 
lensed quasars from tfie statistical properties of tfie unlensed objects in tfie catalog. In §3 we 
summarize tfie tfieory of gravitational lens statistics and discuss tfie effects of tfie selection 
model on tfie probability tfiat a lens is detectable. In §4 we review and expand tfie PS 
lensing model developed by Narayan & Wfiite (1988). In §5 we examine tfie standard CDM 
model, and some of tfie variants suggested to correct tfie problems in COBF normalized 
CDM. In §6 we consider sources of systematic error in tfie calculation and fiow tfiey limit 
tfie cosmological constraints, and in §7 we review tfie results and discuss tfie requirements 
for better wide separation lens surveys. 

2 A Selection Effects Model for Heterogeneous Quasar Catalogs 

We examine tfie statistics of large separation lenses in tfie IIB93 catalog, and Table 1 
summarizes all tfie known lensed pairs, candidate pairs, and associated quasars witfi sepa- 
rations larger tfian 3'.'0 in tfie IIB93 catalog or commonly appearing in lists of gravitational 
lenses. Tfie known galaxy scale lenses witfi separations smaller tfian 3'.'0 (see review by 
Surdej & Soucail 1994) and tfie eigfit objects in Table 1 are tfie full sample of objects in tfie 
IIB93 catalog witfi redsfiifts above 1.0, redsfiift differences smaller tfian 0.01, and separations 
smaller tfian 1'.^. 

Wide separation lenses are always resolved, so given one image of tfie lens witfi magni- 
tude nil brigfiter tfian tfie magnitude limit ?7i;, tfie selection function requires tfiat any second 
lensed image witfi magnitude m2 must also be brigfiter tfian tfie survey magnitude limit. Tfie 
magnitude limit varies from survey to survey, but we would like to fiave a plausible model for 
tfie average dynamic range between a quasar in a fieterogeneous catalog and its magnitude 
limit. We can do tfiis based on two plausible assumptions: (1) all surveys fiave magnitude 
limits, and (2) surveys for fainter quasars do not substantially overlap surveys for brigfiter 



^Be warned, however, that a search of the HB93 catalog with our selection criteria will find seventeen ad- 
ditional pairs. Six of these "pairs" (two pairs near NGC 3384 (Q1045-I-128, Arp et al. 1979), one pair 
near NGC 2683 (Q0849-F336, Arp 1983), and three near Q1549-F486 (Arp & Surdej 1982)) have sepa- 
rations much larger than one arc-minute. The discovery papers did not give separate positions for the 
quasars in the field, so all of the quasars in each field were given a common position. Eleven other 
quasars (Q0057-352, Q0059-411, Q0252-F016, Q0255-015, Q0851-M97=LB8863, Q0953-F549, Q1128-M05, 
Q1151-F068, Q1208-M42, Q1209-M07, and Q1213-M55) appear to be double entries of the same quasar. 



Table 1: Quasar Lenses and Pairs With Separations Larger Than 3'.'0 



Name 


Zs 


A^ 


niA 


niB 


Band 


Lens? 


Why 


Who 


Q0957+561 


1.41 


6'.'1 


17.5 


17.7 


B 


Yes 


Q 


Walsh et al. 1979 


Q2016+112 


3.27 


3'.'6 


22.9 


23.2 


i 


Yes 


L 


Lawrence et al. 1984 


Q1120+019 


1.46 


6'.'5 


16.2 


20.8 


B 


? 


L 


Meylan & Djorgovski 1989 


Q1429-008 


2.08 


5'.'1 


17.7 


20.8 


R 


? 


L 


Hewett et al. 1989 


Q1635+267 


1.96 


3'.'8 


19.2 


20.8 


B 


? 


Q 


Djorgovski & Spinrad 1984 


Q2345+007 


2.15 


7'.'3 


19.5 


21.0 


B 


? 


Q 


Weedman et al. 1982 


PKS1145-071 


1.35 


4'.'2 


18.0 


18.8 


B 


No 


L 


Djorgovski et al. 1987 


Q1343+264 


2.03 


9'.'5 


20.8 


20.9 


B 


No 


Q 


Crampton et al. 1988 



Notes: The entries in the Lens? column are Y if the object 
N if it is generally believed not to be a lens, and ? if its 
in the Why column are Q if the object was found as part 
object was found as part of a lens survey. 



is generally believed to be a lens, 
status is uncertain. The entries 
of a quasar survey, and L if the 



quasars. The first assumption is trivial, but the second requires some justification. The 
surface density of quasars is a steep function of magnitude (the surface densities of quasars 
brighter than 15, 19, and 21 B mags are 1.7 X 10""^, 4.3, and 33 per square degree for the 
redshift range {) < z < 2.2 (Hartwick & Schade 1990), so to find equal numbers of quasars, 
bright quasar surveys cover large areas and faint quasar surveys cover small areas. Since the 
total area surveyed for faint quasars is much smaller than that surveyed for bright quasars, 
the typical bright quasar is not part of a faint quasar survey. Thus the magnitude limit for 
finding companions to any quasar is determined by the magnitude limit of the survey that 
found the quasar. 

We consider a model survey for quasars with a limiting magnitude nii. We assume that 
the survey finds all quasars brighter than the limiting magnitude in a region much larger 
than the largest interesting lens separation. The model neglects photometric errors and 
Eddington bias (see Hartwick & Schade 1990). We model the quasar apparent magnitude 
number counts by a broken power law 



^ _ /V / 10"("-"") m < mo 

dm~ ° 1 lO'^^"-"") m>mo 



:i^ 



where a ~ 1.12, fi ~ 0.18, and niQ ~ 19.1 B mags (Hartwick & Schade 1990, Boyle et 
al. 1990, Wallington & Narayan 1993). The apparent magnitude of the break niQ is nearly 
constant for the redshift range 1 < z < 3.5. In this paper we are uninterested in the absolute 
normalization of the number of quasars A^o- The probability distribution for the magnitude 
difference A?7i = nii — m between a survey quasar with magnitude m and the magnitude 
limit is 



dP 

dAm 



^i dN , 

— — dm 
X) dm 



— (?71; — A?7ll 

dm 



for 



Am > 0. 



For surveys with magnitude limits brighter than the break magnitude, nii < ttiq, the differ- 



ential and integral probability distributions take the simple forms 
dP 



dAm 



[a In 10] 10""^" and P(< Am) = 1 - lO""^". (3) 



In this limit, the mean dynamic range is (Am) = (alnlO)"^ = 0.34 mag, the median is 
Ami/2 = a~^ log 2 = 0.24 mag, and 90% of the quasars have Am < a~^ = 0.79 mag. The 
dynamic range between a quasar and the magnitude limit in any magnitude limited bright 
quasar survey is very small. When the magnitude limit becomes fainter than the break 
magnitude ttiq, the average dynamic range increases. In the limit that nii ^ niQ the mean 
dynamic range is (Am) = (/31nlO)~^ = 1.55 mag, the median is Ami/2 = /3~^log2 = 1.67 
mag, and 90% of the quasars have Am < /3~^ = 5.5 mag. For comparison, surveys for lensed 
quasars (eg. Maoz et al. 1993, Surdej et al. 1993) have typical dynamic ranges of 4 to 5 
magnitudes. 

We test this model by examining the magnitude differences between pairs of quasars in 
the IIB93 catalog. We took the ~ 5000 quasars with measured V magnitudes and 1 < z < A 
in the IIB93 catalog, and found the nearest neighbor for each quasar excluding known lenses 
and the pairs in Table 1. Figure 1 shows the mean and the dispersion of the magnitude 
differences of the pairs as a function of the V magnitude of the brighter quasar. The pairs 
are collected in bins one magnitude wide that are subdivided into a maximum of three 
smaller bins as the number of quasars in the bin increases. Figure 1 also shows the mean 
and dispersion predicted for Am by the selection function model assuming each quasar is 
at the median magnitude for its magnitude limit, and an average B — V color of 0.2 mag. 
Thus, a 18 V mag quasar is assumed to come from a survey with a magnitude limit of 
?7i; = 18 + a~^ log 2 = 18.24 V mag. The agreement between the data and the model 
for both the average magnitude difference and the spread in the magnitude difference is 
remarkably good. 

Figure 1 also shows the locations of the eight quasar pairs in Table 1, labeled by the type 
of survey that found the pair. As expected, the large magnitude difference pairs (Q1120+019 
and Q1429— 008) were found as part of a lens survey. The four pairs found in the quasar 
surveys roughly follow the expected selection function. Although PKS 1145—071 was found 
in a lens survey, it was selected because its image was visibly elongated in the original quasar 
finding chart (Djorgovski et al. 1987). This explains why it lies in the range detectable by 
the original quasar survey. The distribution of lensed pairs will not have the same statistical 
properties as the distribution of unlensed pairs because of magnification bias. 

3 The Effect of the Selection Function on Finding Lenses 

We examine the effect of the selection model on the expected number of lenses using 
the simple singular isothermal sphere (SIS) model for the lenses. A SIS lens with velocity 
dispersion a produces two images with angular separation 9>tt{(t j c^Dls I Dqs where D^s 
and Dos are the proper motion distances between the lens and the source and the observer 
and the source respectively. The integral probability distribution for the two images having 
a total magnification larger than M is P{> M) = A/M'^ with Af > 2, and the fiux ratio 
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Fig. 1. -Average (points) and rms scatter (vertical error bars) of nearest neighbors with redshifts 
l<z<4m the HB93 catalog as a function of magnitude The horizontal error bars show the width 
of the magnitude bins. The lines show the predicted average (solid) and rms (dashed) width of the 
distribution predicted by the model selection function. The quasar pairs in Table 1 are marked by 
a Q if the pair was found in a quasar survey, and an L if the pair was found in a lens survey. 
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Fig. 2. -Magnification bias as a function of B magnitude. The left panel shows the bias factor if 
we find all lenses (heavy solid), all lenses with magnitude differences smaller than the median for 
the selection function (light solid/points), and the average of the bias factor over the distribution 
of magnitude differences (light solid). The right panel shows the completeness, or the ratio of the 
magnification bias with a selection function to the magnification bias when we find all lenses. 



between the two images is / = (M - 2)/(M + 2) with < / < 1 (Gott & Gunn 1974, 
Turner, Ostriker, & Gott 1984). The integral probability lor finding two images with a fiux 
ratio larger than / (closer to unity) is 



Pi>f) 



1-/ 
1 + / 



(4) 



For bright quasars the median dynamic range in our model lor the selection lunction is 
Ami/2 = 0.24, so we only find lenses with f > fi = — 2.51og A?7ii/2 = 0.80 and P{> fi) = 
0.012. II we average P{> f) over the distribution ol dynamic range (eqn. 2), we find an 
average ol P{> fi) = 0.043 ol the lensed bright quasars. For comparison, Wambsganss et 
al. (1994) assumed a uniform dynamic range ol 1.5 magnitudes, lor which P{> fi) = 0.36. 
In the SIS model, this assumption overestimates the expected number ol bright lenses by a 
lactor ol 7. II we relied on the lens cross section to determine the number ol lenses, lensing 
would be useless as a cosmological test. 

Fortunately, predictions ol the number ol lenses found in real surveys must include the 
effects ol magnification bias to correct lor the difference between the number ol quasars at 
the magnitude ol the unlensed source and the number ol quasars at the magnitude ol the 
lensed source (Gott & Gunn 1974, Turner 1980). The bias lactor is 



Bim] 



dN 
dm 



/l'^*'Si5£('" + -"'°«*') f-^) 



where m is the magnitude ol the quasar, M is the total magnification, and Mum is the 
magnification at which the images have the minimum detectable fiux ratio (Fukugita & 
Turner 1991, Kochanek 1991). The completeness, or traction ol the lenses we can detect at 
a given magnitude, is the ratio ol the bias lactor with Mum = 2(1 + fi)/{l — fi) and the 
bias lactor with Mum = 2. Figure 2 shows the magnification bias assuming we find all lenses 
BT^m), all lenses with fiux ratios larger than the median dynamic range predicted by the 
selection lunction, Bii2{m), and the average ol the magnification bias over the probability 
distribution lor the dynamic range (eqn. 2), {B)[m). The magnification bias compensates lor 
the low optical depth at most magnitudes, although the Cen et al. (1994) model significantly 
overestimates the amount ol bias. The standard quasar number counts used by Cen et al. 
(1994) based on Fukugita & Turner (1991) also have a significantly shallower slope lor the 
number counts ol the bright quasars [a = 0.86, /3 = 0.28) than more recent determinations 
(a = 1.12, /3 = 0.18; see Boyle et al. 1990, Wallington & Narayan 1993). 

Given the selection lunction lor finding lensed quasars and the magnitude limit aver- 
aged magnification bias lactor {B)[m)^ the probability that a quasar with redshilt z and B 
magnitude m is a lens with an image separation larger than A^ is 



Dos Dl^dDoL fDLsV 



CO 



Jo (1 + ilKP>^j^rfj y'^ \P>osy Jcrmu \cJ da 

, , , (^) 

where dn/dcr[cr^ z) is the comoving number density ol lenses with velocity dispersion a at 

redshilt z^ Umin = c{A9Dos/8TrDLsY is the smallest velocity dispersion that can produce 

7 



an image separation A^ at redshift z^ and c is the speed of light. The distances Dqs^ Dql^ 
and Dls are the proper motion distance to the source, to the lens, and between the lens and 
the source. The Hubble radius is r^ = c/i^o, Hq = lOOh km s~^ Mpc~^, and ftx = 1 — Oq — Aq 
is the "curvature density" for a cosmological model with matter density Oq and cosmological 
constant Aq (Carroll et al. 1990, Kochanek 1993a). The distance D^s can be determined 
from DoL and Dqs by the relation Dls = Dos{l + ^kDIlVh^Y'^ - Dol{1 + ^RDlsrH^f'^ 
From equation (6) we can also compute the probability that a lens has separation A^ to be 

dP vr^^,„,, , [Dos Dl^dDoL dn 



where a = c{A9Dos /SttDlsY ■ The integrand of dP/dA9 (converted to a differential with 
respect to lens redshift) is the lens redshift probability distribution for lenses with image 
separation A^. 

Given the number density of potentials dn/da we compute the probability pi that the 
i^^ quasar with B magnitude rrii and redshift Zi in the HB93 catalog is lensed and detectable 
given our model selection function. We use the maximum likelihood formalism of Kochanek 
(1993b) to estimate the likelihood that different models fit the lens data. For a sample of 
Nu unlensed quasars and Nl lensed quasars the likelihood L of the observations is 

lnL = -J2Pk + T.^nP^ (8) 

fc=l 8=1 

where we use the expansion ln(l — Pk) = —Pk because pk <^ I. If L^ax is the maximum 
value of the likelihood for some range of model parameters, then the function —2 ln[L / L^ax) 
is asymptotically distributed like the x^ distribution (Lupton 1993). The 68% (Icr), 90%, 
95.4% (2(t) and 99% confidence levels on parameters in one (two) dimensions are where the 
likelihood is 60.7% (31.7%), 25.8% (9.98%), 13.5% (4.57%), and 3.63% (1.00%) of the peak 
likelihood. 

We present two standard calculations. The first is the likelihood of finding lenses larger 
than A^ = 5'.'0, and the second is the likelihood of finding lenses larger than 3'.'0 including 
the likelihood of finding the lenses with their observed separations (eqn. 7 instead of eqn. 
6). By examining these two different statistical estimates we check whether the conclusions 
are sensitive to the lower angular cutoff (5'.'0 versus 3'.'0) and whether the results are strongly 
sensitive to the observed separations. We know that Nl is at least equal to one, because 
there is one certain lens in the sample, Q0957+561. There might be as many as three lenses 
larger than 5'.'0 if we include Q2345+007 and Q1343+007. We generally present results for 
Nl = 1 and Nl = 3 to show the sensitivity of the conclusions to the ambiguities in the 
numbers of lenses. When we use the 3'.'0 angular limit we drop the Q1343+007 pair and 
replace it with the lens candidate Q1635+267. 

4 Estimates of the Number Density of Potentials 

Following the approach of Narayan & White (1988) we use the PS model to estimate 
the comoving number density of clusters dn/da as a function of redshift z and velocity 

8 



dispersion a. Let 6c{z) be the critical overdensity Sp/{p) that can collapse before red- 
shift z^ and assume that the fluctuations at a comoving scale of Tq are Gaussian with an 
rms linear overdensity of A(ro). The fraction of the universe in such high density regions 
is F{ro,z) = {2TrA{ro))~^''^ f^r^exp^—u"^ /2A{roY)du. PS suggested that the number of 
clumps with initial comoving radii in the range Tq to Tq + dro can be approximated by 
f[ro^ z)dro = —2[dF/dro)dro^ where the 2 is inserted to ensure that J^ fdro = 1. Pertur- 
bations exceeding the critical overdensity collapse and virialize. Bond et al. (f99f) give a 
rigorous derivation of the PS ansatz. After multiplying /(tq, z) by the appropriate Jacobian 
and normalizing it, we find that the comoving number density of halos at redshift z is 



dn 
da 



cr,z 



-3Sc{z) din A 
(27r)3/2r3aAJlnro 



exp 



^li^) 



2A2 



(9) 



if the velocity dispersion of the collapsed object is cr oc Tq. If the power spectrum of the 
fiuctuations is \6k\'^ then the variance of the fiuctuations on scale Tq is given by the convolution 



A\ro) 



f27rl 



47rP dk\Sk\^W^kro) 



where 



W{x) 



smx 



X cos X 



(10) 



is the Fourier transform of the top-hat window function. We normalize the power spectrum 
by the rms fiuctuation erg = A(ro = 8/i~^Mpc) on the scale of Tq = 8h~^ Mpc, or by the 
"bias" factor b = I /erg. For any cosmological model and power spectrum we can compute 
the expected number of lenses given the function 6c{z) for the smallest perturbation that 
can collapse at redshift z and the relationship between a and Tq. 

Bartelmann et al. (1993) and Lacey & Cole (1993) show that the redshift z at which an 
initial perturbation with amplitude Si at redshift Zi ^ Oq ^ collapses in a Aq = cosmological 
model is 



1 



u^y^ 


' 5Oo(l + ^0 


where the age of the universe at redshift z is 


T{z) = H-'F,{z) = H-' 


(l + Oo^)i/2 
[(1-Oo)(l + ^) 



ttOo 



2/3 



+ I-O0 



111 



Oo 



2(1 - Oo)3/2 



cosi 



2(1 -Op) 
Oo(l + ^) 



+ 1 



12) 



We use a power spectrum that is normalized by the fiuctuations today, so we must use linear 
theory to convert from bi to b^. The growing mode changes its amplitude with redshift by 
the function (Peebles 1980) 



DM 



1 + - + 

X 



3 3(l + x)i/2 



p3/2 



ln[(l + 



Nl/2 



.1/21 



where 



Oo' 



1 



1 + ^ 



so 



^ ^ 2 1 - Oo 



ttOo 



2/3 



+ 1-O0 



(13) 



(14) 



In the limit that Oq ^ 1 the function i)i[Oo^ - 1]/(1 - Oq) -^ (2/5) and Fx{z) -^ (2/3)(l + 
z)~'^''^, and we recover the standard result that Sc{z) = (3/5)(37r/2)2'"3(l -|- z). 
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The mass of the collapsing perturbation is Af = ^tt por^j?) where po = 3HQfto/8TrG is 
the average mass density. The object collapses and virializes to form an isothermal sphere 
of velocity dispersion cr, radius i?j, and a mass inside Rj of Af = 2cr^Rf /G. Equating the 
two masses we find that 

1/2 

(15) 



a = Horo9.l'^ 



To 



4R 



The virial theorem for the collapsing perturbation shows that the final virialized radius is 
Rf = Rmaxl'^ (Lahav et al. 1991), where Rmax = ?"o(l + Zi)~^{b8il'i — ei)~^ . Combining the 
expression for Rmax with equation (15), we find that 

a = Horonl^'2-'^\'^'Friz)-'^\ (16) 

In the limit that Oq = 1 we find that standard result that a = 2-^I^{?,tt 1 2f I ^Hor q{1 + zfl^. 

For fiat universes with a positive cosmological constant (Oq + Aq = 1) none of the 
integrals needed for the PS model can be done analytically, but following Richstone, Loeb, 
& Turner (1992) we numerically solve the implicit equations for Sc- The age of the universe 
T = Hq Fi{z) is given by 

Fi{z)= f ' u^/2Ju [u^Ao + Oo]"^^^ (17) 

and the amplitude of the growing mode is proportional to 

r . xQ ii/3 /■(i+^)~^ ,/o r Q 1-3/2 , , 

D,{z)=[Xo + {l + zfno\' J u3/2ju [u^Ao + Oo] . (18) 

The collapse time Tc{6i) for a perturbation with fractional overdensity Si at redshift Zi ^ 1 
is 

HoT,{6,) = 4Ao ^/^ r sin^ 6 [u2/ui - sin^ o] ~^'^ [sin^ Q - u^ju^ ~^'^ (19) 

where M3 < < Ui < M2 are the roots of the cubic polynomial in u = r/ro 

= 1 - e, + (5, - 5(5,m/3 + e^M^ (20) 

and ti = Ao [Aq + Oj(l + ZiY] ^ 1 is the cosmological constant or the deviation of the matter 
density from unity (Oj- = 1 — €{) at redshift Zi. The root Ui = Rmax/i'o is the ratio of the 
maximum radius Rmax of a perturbation to its initial radius Tq, and the perturbation must 
be larger than Si ^ (9/5)(ej/4)^'"^ to collapse in finite time. The function Sc{z) is constructed 
by solving the implicit equation Fi[z) = HoTc{6i) for 6i[z) and using the growth rate of the 
growing mode to convert from Si to 6c{z) = [Di[0) / Di[zi)] 6i[z). 

The velocity dispersion of the collapsed object in the models with a cosmological con- 
stant must correct for the change in the virial radius from the energy associated with the 
cosmological constant (Lahav et al. 1991). If ?y = Xj^irGpta is the ratio of the cosmological 
constant to the average density when the perturbation turns around and collapses, then 
the ratio of the virial radius Rj to the turn around radius Rmax satisfies the cubic equa- 
tion 2r][Rf / RmaxY — (2 -|- r])[Rf / Rmax) + 1=0. When rj = we find the standard result 
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Rf I Rmax = 1/2, but when Aq > 1 the lack of the repulsive force from a positive cosmological 
constant in the collapsed halo leads to a smaller virialized object. The solution for the virial 
radius is well approximated by Rf j Rmax = 0.5 — 0.138?7 + 0.0034?y^ for ?y > 0. (For ?y > 
this is more accurate than the approximation in Lahav et al. (1991).) 

5 Lensing in Different Cosmological Scenarios 
We use the approximate power spectrum (eg. Efstathiou, Bond, & White 1992) 

where a = {6A/T)h-^ Mpc, b = {3.0/T)h-^ Mpc, c = {1.7/T)h-^ Mpc and ly = 1.13. 
The fitting formula is valid if the baryon density is much smaller than the cold dark matter 
density. The power spectrum determines the rms fiuctuations A on scale Tq through equation 
(10). We always normalize the power spectrum to the fiuctuations erg on Tq = 8h~^ Mpc, so 
A(ro) = cr8A(ro) and A(8/i~^Mpc) = 1. The calculation depends on the Hubble constant /i, 
the power spectrum normalization CTs, the shape of the power spectrum (n and F), and the 
cosmological model. Before considering individual models, we consider the sensitivity of the 
lens calculation to the various parameters. 

The lens calculations are insensitive to the value of the Hubble constant. The number 
density of galaxies as a function of velocity dispersion (eqn. 9) depends on h^^ but it is 
multiplied by the h~^ dependence of the comoving volume element in the probability calcu- 
lation (eqn. 6). The Hubble constant dependence of the a, 6, and c coefficients of the power 
spectrum, excluding the parameter F, gives the integral a dependence on l/[hro) oc 1/cr, 
again canceling the dependence on the Hubble constant. When we normalize the spectrum 
by (Tg we remove the Hubble constant dependence of the normalization constant B. Thus 
the expected number of lenses depends on the Hubble constant only through changes in the 
shape parameter F, and this dependence is weak for 0.5 < h < 1.0. The COBE estimates for 
(Tg are sensitive to the Hubble constant (eg. Efstathiou, Bond, & White 1992, Bond 1994), 
so we generally show the lens estimates for h = 0.5 compared to the COBE estimates for 
h = 0.5 and h = 1.0. 

When we study cluster lenses with separations larger than 5'.'0, we are examining objects 
with velocity dispersions between 300 km s~^ and the largest objects that can collapse and 
virialize before the present epoch. The velocity dispersion a for a perturbation of scale Tq is 



a{ro) = 950 



ro 



8/i-iMpc 



'l + zY^'kms-' (22) 



in an Oo = 1 cosmology (eqn. 16). Thus the scale used to normalize the power spectrum 
(ro = 8h~^ Mpc) is comparable to the scale producing the lenses, and the number of large 
separation lenses is controlled by dn/da oc ag^ exp{—S^/2al) at the 8h~^ Mpc scale. The 
expected number of lenses depends exponentially on the normalization when erg < 8c, and the 
number varies as a power law when CTs > 6c. Since the critical overdensity is (Tc = 1.69(1 + z) 
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for Oo = 1, the number of lenses will vary exponentially with erg. The smaller separation 
lenses produced by galaxies correspond to smaller, more nonlinear perturbations where dn/da 
is in the power-law regime. 

For a fixed value of erg the number of lenses depends weakly on the shape of the power 
spectrum. The slope of the power spectrum controls the the separation distribution of the 
lenses. Power spectra that increase more rapidly for smaller wavenumbers will produce 
lens separation distributions that decline more steeply with increasing image separation 
than fiatter power spectra. The maximum likelihood method can be made sensitive to the 
separation distribution by using the probability that a lens has a given separation (eqn. 7) 
in the lensed term of the likelihood (eqn. 8) instead of the probability that it is a lens (eqn. 
6). Given the small number of lenses in the current sample this will not strongly constrain 
the shape of the power spectrum, but we include the separation probability distribution in 
one of two standard statistical models to get a feeling for how strongly it will constrain the 
slope of the power spectrum. 

5.1 Standard CDM 

In standard cold dark matter (CDM) models, the primordial scale-invariant spectrum 
is \6k\'^ = Bk [n = 1), the cosmology is fiat with Oq = 1, and the parameter T = h (eg. 
Efstathiou et al. 1992). The only free parameters are the normalization constant B which 
we parametrize by the variance in the density erg on the scale of Tq = 8h~^ Mpc and the 
Hubble constant h. The measurements of the microwave background fiuctuations by the 
COBE DMR experiment (Smoot et al. 1992, Wright et al. 1994, Gorski et al. 1994) require 
-PS = (19.9 ± 1.6)/iK if n = 1 which corresponds to os — (2.8 ± 0.2)h. 



' rms- 



Wambsganss et al. (1994) numerically calculated the optical depth to lensing for a 
standard CDM model with erg = 1.05 and h = 0.5. They found that the optical depth 
for lenses larger than lO'.'O (5'.'0) with fiux ratios smaller than Am = 1.5 magnitudes are 
0.0007 (0.0008), 0.0014 (0.0019), and 0.0020 (0.0027) for source redshifts of z^ = 1, 2, and 
3 respectively. The PS calculation with the same parameters and selection function yields 
cross sections of 0.0004 (0.0006), 0.0012 (0.0017), and 0.0019 (0.0027). The differences in 
the optical depths are 43% (31%), 21% (11%), and 10% (0%). Narayan & White (1988) 
used an approximation for A(ro) that accounts for most of the differences noted by Cen et 
al. (1994) in their comparison to the PS model. Cen et al. (1994) predict 46 lenses larger 
than 4'.'0 in the older Hewett & Burbidge (1989) catalog for erg = 1 including magnification 
bias, comparable to the 27 predicted here even though the samples, quasar number counts, 
and selection model are different. 

The separation distribution of lenses (see Figure 9) in the PS simulations declines more 
slowly than the numerical simulations (for erg = 1.05), but the distributions are similar for 
A^ <^ 40". Detailed comparisons are difficult because Wambsganss et al. (1994) use the 
largest detected image separation, and this procedure will transfer power from large sepa- 
rations to smaller separations because of the five image systems with only two detectable 
images (see §6). The separation scale of the deviations roughly corresponds to perturbations 
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Fig. 3. -Lenses in standard CDM. The left panel shows the expected numbers of lenses larger 
than S'.'O, lO'.'O, IS'.'O, and 20'.'0 as a function of a^ for h = 0.5. Standard CDM normalized to 
COBE requires ag = (2.8 ± 0.2)h. The right panel shows the likelihood ratio as a function of ag for 
lenses larger than S'.'O. The heavy solid curve assumes Q0957+561 is the only lens, the light solid 
curve assumes Q0957+561 and Q2345+007 are lenses, and the dashed curve assumes Q0957+561, 
Q2345+007, and Q1343+264 are lenses. The horizontal line shows the 2a limit on the likelihood 
ratio, and the point labeled COBE shows the best COBE normalization for ag (h = 0.5) and its 
1(7 error bar. 
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Fig. 4. -Lens redshifts in standard CDM. The two panels show the probability distributions for 
the lens redshift given the observed image separations in Q0957+561 and Q2345+007. The dashed 
lines use the COBE normalized values of ag and the solid lines have ag = 0.5. The lines with 
no points use a Hubble constant of h = 0.5, and the lines with points use a Hubble constant of 
h = 1.0. The vertical lines mark the known lens redshift in Q0957+561 and a suggested redshift 
for Q2345+007 (Fischer et al. 1994). 
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with comoving scales significantly larger than the size of the high resolution numerical simu- 
lations [5h~^ Mpc). The integral redshift distribution of the lenses in the PS model declines 
more rapidly at low redshifts than in the numerical simulations. This can be explained by 
adding a core radius comparable to the grid resolution of the numerical simulations [10h~^ 
kpc) in our lens potential model. We consider the effects of ellipticity and core radii in more 
detail in §6. Nonetheless, the two calculations are everywhere in agreement by a factor of 
two or better when there is an appreciable cross section for lensing and on scales smaller 
than 40'.'0. Moreover, the exponential dependence of the number of lenses on erg means that 
large uncertainties in the optical depth, the selection function, and the number of lenses lead 
to much smaller uncertainties in the value of erg. For the standard CDM model, changing 
the expected number of lenses by a factor of two near erg = 1.4 changes the estimate of erg 
by 0.4, and changing it by a factor of two near erg = 0.5 changes the estimate of erg by less 
than 0.1. We discuss the sources of systematic errors in §6. 

Figure 3 shows the expected number of lenses as a function of erg in standard CDM. 
The number of lenses rises exponentially with the bias I/cts, and 60 lenses larger than 5'.'0 
are expected at the COBF normalization of erg = 1.4 it 0.1 for h = 0.5. This is coincidentally 
similar to the cross section calculations of Wambsganss et al. (1994) because an order of 
magnitude error in the estimated cross section is balanced by an order of magnitude error 
from neglecting the magnification bias. Figure 3 also shows the likelihood ratio as a function 
of erg if we include the one real lens (Q0957-|-561), the real lens and the best candidate 
(Q0957-|-561 and Q2345-I-007), or the real lens, the best candidate, and the rejected candidate 
(Q0957+561, Q2345+007, and Q1343+264). The 2(7 range for ag is 0.27 < ag < 0.63, setting 
the lower limit using one lens and the upper limit with all three. The maximum likelihood 
values are erg = 0.40, 0.46, and 0.50 for one, two, and three lenses. These limits are little 
changed by reducing the separation limit to 3'.'0 or changing the Hubble constant to h = 1. 

The lens results rule out standard CDM normalized by COBF. The best fit value for 
erg is consistent with other estimates based on number density of clusters (Peebles et al. 
1989, Frenk et al. 1990, Bahcall & Cen 1992, 1993) or the correlation function of galaxies 
on large scales (Maddox et al. 1990, Picard 1991, Vogeley et al. 1992, Loveday et al. 1992). 
It is, however, an independent test controlled by the properties of the cluster distribution 
at moderate redshifts rather than today. Figure 4 shows the lens redshift probability dis- 
tributions for Q0957-|-561 and Q2345-I-007. The lens redshift in Q0957-I-561 is known to 
be zi = 0.36, and it lies at the peak probability of the probability distribution normalized 
to produce the observed number of lenses [erg ~ 0.5). The COBF normalized distributions 
extend to much higher redshifts because with larger values of erg the clusters form earlier. 
In Q2345-|-007, Fischer et al. (1994) suggest that the lens redshift is zi = 1.48 because there 
are many metal absorption features at that redshift. In COBF normalized CDM zi = 1.48 
is a plausible cluster redshift, but in the models that predict the observed numbers of lenses 
it is extremely unlikely. 
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5.2 Tilted CDM 

The standard scale-free inflationary model assumes that the primordial power spectrum 
is \6k\'^ oc k. One set of models designed to solve the problems standard CDM has in 
simultaneously fltting the COBE scale and smaller scale structures, "tilt" the CDM spectrum 
by making the primordial spectrum \6k\'^ oc k^ with n ^ 1 (eg. Vittorio et al. 1988, Cen et 
al. 1992). Current flts to the COBE data (Wright et al. 1994, Gorski et al. 1994) flnd that 
n = 1.2 lb 0.3, consistent with n = 1. The fltted value for the Q^ms-ps parameter varies with 
the value of n, and Gorski et al. (1994) give Q^ms-ps = 18.7(1 + 0.05n) exp(0.73(l — n)) //K 
for the best flt value as a function of n. Combined with the fltting formula for estimating erg 
from Qrms-ps by Bond (1994) we flnd 

(78 = (2.8 ± 0.2)(r - 0.03)(1 + 0.05n) exp [-1.9(1 - n)] . (23) 

This normalization neglects any tensor-mode contributions from gravitational waves. Eigure 
5 shows contours of the likelihood of fltting the lens data as a function of erg and the exponent 
of the tilt n. The primordial spectrum must be flatter than n = 1 for the lens estimate of erg 
to agree with COBE, simply because a flatter spectrum gives less small scale power on the 
(Tg scale when the larger scale is flxed by COBE. 

Eor h = 0.5 the lens and COBE limits agree if 0.3 ^ n ^ 0.7 (2(t), and if h = 1.0 they 
agree if 0.0 ^ n <^ 0.3 (2(t). The lens calculations were done with h = 0.5, but they are 
insensitive to the difference between h = 0.5 and h = 1.0, so we can effectively make the 
comparison using the h = 0.5 calculation. The tilt exponents at the upper end of the allowed 
range for h = 0.5 are marginally consistent (at 2cr) with the COBE limits on the value of n, 
although adding a tensor-mode/gravitational wave contribution to the COBE signal would 
improve the agreement. The left panel of Eigure 5 only includes lenses with separations 
larger than 5'.'0. Since the likelihood has no information other than the expected number 
of lenses, it cannot differentiate between different values for the exponent n without outside 
information on erg. Tilting the spectrum rearranges the distribution of image separations, 
and flat spectra (n < 1) produce relatively more large separation lenses than steeper spectra 
(n = 1). If we normalize erg so that we flnd three lenses larger than 5'.'0, then for n = 0, 
n = 1, and n = 2 we expect 0.6, 0.3, and 0.1 lenses larger than lO'.'O (see Eigure 9). 

We can add sensitivity to the value of n by using the probability that the lenses have 
their observed separations, and the right panel of Eigure 5 uses the modifled likelihood for 
the lenses larger than 3'.'0. By comparing the two estimates we draw three conclusions. Eirst, 
the estimated range for erg does not markedly change when we alter the angular selection 
function. Second, if there is only one observed lens we have no leverage to determine the 
value of n because even for the whole range of < n < 2 the separation of Q0957-|-561 never 
becomes unlikely. Third, if we add even a small number of additional lenses, the likelihood 
contours begin to close. If Q0957-|-561, Q2345-I-007, and Q1635-I-267 are all lenses, then 
values of n smaller than n = 0.4 are ruled out at the la level. When n ~ the expected 
separation distribution is too flat for all three lenses to have separations between 3'.'0 and 
7'.'0 (see Eigure 9). 
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Fig. 5. -Lenses in tilted CDM. The left panel uses the likelihood of finding lenses with separations 
larger than S'.'O. The solid lines are likelihood contours for finding one lens (Q0957+561) and 
the dashed lines are likelihood contours for finding three lenses (Q0957+561, Q2345+007, and 
Q 1343+264) as a function of the exponent n and the ag normalization. The right panel uses the 
likelihood of finding lenses with separations larger than 3'.'0 and the likelihood that the observed 
lenses have their measured separations. The solid lines are likelihood contours for finding one lens 
(Q0957+561) and the dashed lines are likelihood contours for finding three lenses (Q0957+561, 
Q2345+007, and Q 1635+267) as a function of the exponent n and the ag normalization. The 
contours are the la and 2a confidence intervals of the likelihood ratio for one parameter. The 
best fit models are bounded by the contours. The heavy solid lines show the COBE normalized 
estimates for ag when h = 0.5 and h = 1.0. 
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Fig. 6. -Lenses in fitted CDM. The left panel uses the likelihood of finding lenses with separations 
larger than 5'.'0, and the right panel uses the likelihood of finding lenses with separations larger 
than 3'.'0 and the likelihood the observed lenses have their measured separations. The heavy solid 
lines show the COBE normalized estimates for ag, and the cases and contour levels are the same 
as in Figure 5. 16 



5.3 The Hubble Constant And Varying T 

In the standard power spectrum with fto = 1^ T = h depends only on the Hubble 
constant. This is the only place where the calculation depends on the value of the Hubble 
constant. We extend the calculations over a wider range for F than is plausible for the 
Hubble constant, and in this broader treatment F is just a fitting parameter in the power 
spectrum. Figure 6 shows the dependence of the likelihood on F and erg for both A^ > 5'.'0 
with the total lens probabilities and A^ > 3'.'0 with the lens separation probability. When F 
is increased from 0.5 to 1.0 the best fit values of erg shift downwards by 0.1. Since the shift 
in (Tg over the range 0.5 < h < 1.0 for the Hubble constant is significantly smaller than the 
error bars, we can regard lens models with h = 0.5 as being equivalent to the same model 
with h = 1.0. There is no plausible value of the Hubble constant that can bring the lens 
and COBE estimates into agreement, but if we regard F simply as a fitting parameter, they 
are consistent when 0.15 ^ F ^ 0.30 (2(t) and 0.3 ^ (Jg ^ 0.7. These ranges are similar 
to the ranges found from the correlation function on scales near 10h~^ Mpc (Maddox et al. 
1990, Picard 1991, Vogeley et al. 1992, Foveday et al. 1992) and from observed numbers of 
large clusters (Bahcall & Cen 1992, 1993). The lens estimate for erg as a function of F is a 
measure of the rms perturbations in the mass distribution, so the agreement with estimates 
of (Ts from luminous objects implies that there is little or no bias on cluster scales. With the 
inclusion of the lens separations, the likelihood contours begin to pinch off for small values 
of F (see Figure 9). 

5.4 Open and Flat Low Oq Universes 

Reducing the matter density has three different effects on lens statistics. The first effect 
is that lower density universes have larger comoving volumes to a given redshift, so for a 
constant comoving density of objects we expect more lenses. An empty universe has twice 
as many lenses (Turner, Ostriker, & Gott 1984) and a fiat universe with Aq = 1 has ten 
times as many lenses (Turner 1990) as a fiat universe with Oq = 1- This allows us to use the 
incidence of galaxy scale lenses to determine the cosmological model (eg. Kochanek 1993b, 
Maoz & Rix 1993). The second effect is that perturbations stop growing at low redshifts (at 
z ~ Oq for open universes) so clusters must collapse earlier in low matter density universes 
if they are to be seen today. Richstone et al. (1992), Facey & Cole (1993), Bartelmann 
et al. (1993) use this to argue that high values of Oq are needed to explain the observed 
substructure and rapid evolution of clusters at low redshifts. Both of these effects work to 
increase the expected number of lenses. The third effect is that a perturbation on scale Tq 
contains less mass in a low Oq universe, and when it collapses it tends to have a larger virial 
radius. The mass of the perturbation scales with OoTq, and the virial radius is proportional 
to the maximum expansion radius. 

For a fixed value of CTs, the third effect is the dominant one, and lowering Oq reduces the 
number of observable lenses. Because the velocity dispersion of the collapsed perturbation 
for a fixed value of Tq is smaller, lenses with a fixed image separation are produced by 
perturbations on larger and larger comoving scales as Oq is reduced. At fixed CTs, this means 
that the rms fiuctuation amplitude of the perturbations producing a fixed image separation 
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decreases. This makes changes in Oq for a fixed erg similar to changes in erg for a fixed Oq, 
because of the way it changes the power spectrum normahzation on the scales probed by 
lenses with a fixed separation. Since the number of lenses is exponentially sensitive to the 
normalization, the third effect dominates over the other two. This is very different from 
galaxy scale lenses, where the geometric effects of changing the cosmological model are more 
important than evolution (Mao f99f , Mao & Kochanek f994, Rix et al. f994). 

Figure 7 shows the likelihood contours for open universes compared to the COBE nor- 
malizations determined for Oq < f CDM models by Kamionkowski & Spergel (f994) renor- 
malized to fit the Gorski et al. (f994) rms quadrupole estimates. The approximate 2cr limit 
is 0.f5 ^ ftoh <^ 0.3, although this overestimates the upper limit for h = 0.5 by about O.f 
in Oo- Note that when we add the separation probabilities and three lenses, the likelihood 
contours have pinched off at both low Oq giving a nominal one standard deviation limit of 
Oo ^ 0.25. Obviously this is a weak conclusion, but it shows how the separation distribution 
of a larger sample of large separation lenses can constrain the slope of the power spectrum. 
Figure 8 shows the likelihood for fiat universes with a cosmological constant. The limits are 
consistent with the COBF Qrms-ps if 0.45 ^ Aq ^ 0.85 (2(7) for h = 0.5 and 0.85 ^ Aq ^ f .0 
(2(t) for h = f .0. We again see signs that the observed separation distribution may be too 
steep for low values of Oq- Gen et al. (f994) also find that models with Aq = 0.3 have nearly 
an order of magnitude fewer lenses than models with Aq = 0.0 when erg = I. Recall, how- 
ever, that gravitational lensing by galaxies strongly rules out large values of the cosmological 
constant (the two standard deviation upper limit is Aq < 0.6, Kochanek f994), and that the 
observed, low value of the microwave background quadrupole compared to Q^ms-ps supports 
this conclusion (Sugiyama & Silk f994). 

6 Systematic Errors 

Before concluding we should discuss some of the systematic errors in this calculation. 
We can divide the problems into shortcomings with the PS method as a means of estimating 
the number of halos, shortcomings in converting a collapsed perturbation on scale Tq into 
a virialized object with velocity dispersion cr, and shortcomings in modeling the lenses as 
circular, singular isothermal spheres. To illustrate the effects of these systematic errors, we 
examined the standard CDM model to see how the estimate of erg depends on changes in 
each of these parts of the model. 

The critical part of the PS model is Sc{z)^ the critical overdensity collapsing at redshift z. 
While errors in the details of the lens model and the selection function cause only logarithmic 
changes in estimates of CTs, changes in 6c lead to proportionate changes in erg because the 
expected number of lenses depends roughly on exp(— (5^/2(t|). If we vary 6c between 0.5 and 
1.5 of its standard value, then the estimated value of erg changes by ±0.2 for finding one 
lens and by ±0.25 for finding three lenses. Reducing 6c lowers the perturbation amplitude 
needed to produce the observed number of lenses. This scaling is roughly as expected, since 
a 50% variation in 6c produces a 50% variation in the best fit value of erg. If we multiply 
6c by 1.5 and assume a sample with three lenses, then the the 2cr upper limit on erg rises to 
0.92. 
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Fig. 7. -Lenses in low fio CDM. The left panel uses the likelihood of finding lenses with separations 
larger than S'.'O, and the right panel uses the likelihood of finding lenses with separations larger 
than 3'.'0 and the likelihood the observed lenses have their measured separations. The heavy solid 
lines show the COBE normalized estimates for ag from Kamionkowski & Spergel (1994), and the 
cases and contour levels are the same as in Figure 5. 
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Fig. 8. -Lenses in fio + •^o = 1 CDM. The left panel uses the likelihood of finding lenses with 
separations larger than S'.'O, and the right panel uses the likelihood of finding lenses with separations 
larger than 3'.'0 and the likelihood the observed lenses have their measured separations. The heavy 
solid lines show the COBE normalized estimates for ag, and the cases and contour levels are the 
same as in Figure 5. 
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The virial theorem estimates for converting the perturbation scale Tq into a velocity 
dispersion are also crude. The estimate of erg is sensitive to the estimated ratio of the final 
virialized radius because by changing the relation between Tq and a we shift the comoving 
scale that produces lenses of a given separation. Lowering the final virial radius means 
that a given velocity dispersion can be produced by smaller perturbations, reducing the 
estimate of erg. We used the estimate (for Aq = 0) that the final virialized radius was one 
half the maximum expansion radius Rj = Rmax/'^^ and the velocity dispersion estimate is 
proportional to the square root of the ratio a oc (Rmax/ RfY ■ For example, the velocity 
dispersion estimate rises by 10% if we treat the collapsed perturbations as Jaffe (1983) models 
and normalize the velocity dispersion by conserving energy instead of mass. If we halve the 
virial radius Rj = Rmaxj ^ then the velocity dispersion associated with a perturbation is 40% 
larger, and the best fit values of (Tg, are reduced by 0.13 for one lens and by 0.19 for two 
lenses, while if we double the virial radius Rj = Rmax then the velocity dispersion is 40% 
lower, and the best fit values of erg are increased by 0.23 for one lens and by 0.32 for three 
lenses. If we take three lenses and increase the virial radius by a factor of two, the 2cr upper 
limit on erg is 1.14. The value of erg is roughly scaling as [2Rf / R^ax)- 

If we add ellipticity to the lens model we do not change the average cross section of the 
lenses but we radically alter the range of image morphologies because the lenses produce four 
image as well as two image systems. This can alter the number of observable lenses through 
changes in the magnification bias and the detectability of the new image configurations. We 
generated Monte Carlo samples of 10"* lenses produced by singular isothermal spheres in an 
external shear field with dimensionless ellipticity 7 = 0.1 (see Kochanek 1991) at 15, 16, 
17, 18, and 19 B magnitudes. The average magnification bias before applying the selection 
function is nearly identical to the circular model. We then applied the selection function to 
the Monte Carlo catalogs, using the average estimate for the limiting magnitude [{Am)) and 
found that the expected number of lenses was within 20% of the estimates for the circular 
lens. This has a negligible effect on estimates of erg. The introduction of elliptical lenses 
does radically alter the expected morphologies because the four image systems are almost 
always detected - the incompleteness comes entirely from reducing the numbers of two image 
systems. Thus at 16, 17, 18, and 19 B mags the expected fraction of four image lenses is 
87%, 45%, 33%, and 30% of the total, even though they consist of only 63%, 41%, 13%, 
and 4% of the lenses at those magnitudes. The preponderance of the four image systems 
is not too overwhelming to be inconsistent with finding only two images systems in the 
existing sample. Moreover, particularly for the brighter systems, roughly 85% of the four 
image systems are only detectable as image pairs in a quasar survey, with the remaining 
two images lying beyond the magnitude limit. Surveys tend to find the quads as merging 
pairs on a critical line because they are substantially magnified and have similar fiuxes. The 
introduction of ellipticity has no interesting effect on estimates of erg. 

The last systematic problem we consider is the addition of a core radius to the lens 
model. We know from the models of giant arcs (see the review by Soucail & Mellier 1994) 
that the core radii of the arc producing clusters are compact enough to produce well separated 
tangential and radial critical lines for comparatively low redshift sources [zs ^ 1), and we 
know from studies of the galaxy scale lenses that the core radii of galaxies are very compact 
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(eg. Wallington & Narayan 1993, Kassiola & Kovner 1993), but we have no direct information 
on the core radii of groups and clusters with intermediate masses. We emphasize comparisons 
to limits from lensing because they are the most direct constraints on the core radii used in a 
lensing calculation. Numerical simulations of structure formation uniformly lead to collapsed 
objects with core radii set by the smoothing length of the simulation, consistent with all 
objects having nearly singular dark matter halos (eg. Dubinski & Carlberg 1991). When 
we consider adding core radii to the models we must consistently include the modifications 
of the magnification bias, because correcting the expected number of lenses using only the 
changes in the lens cross section produced by core radii leads to serious quantitative errors. 
If the lenses have the density distribution p oc (r^ + s^)~^ then the lenses are subcritical if 

where the D^ are angular diameter distances, s = 10sioh~^ kpc, and b = 4:Tr{a/cYDf^g/ Dq^ 
is the critical radius of the singular model (Hinshaw & Krauss 1987). For Oq = 1 and a 
source at redshift Zs = 2, the core radius must be smaller than s < 28h~^cr^ kpc for the 
lens to produce multiple images. The 10h~^ kpc comoving grid scale in the Wambsganss et 
al. (1994) simulations is dangerously close to this limit. If we assume that the core radii 
of clusters scale with s oc cr^ so that the ratio s/b is independent of velocity dispersion and 
all clusters have an equal capacity for lensing, then for a source at Zs = 2 we find that the 
expected number of lenses is reduced by a factor of two for a core radius of s = 15h~^cr^ 
kpc, it is reduced by a factor of ten for a core radius of s = 19h~^crl kpc, and it is zero 
for s > 28h~^al kpc, including the effects of magnification bias. For intermediate core radii 
of order 7h~^cr^ kpc the expected number of lenses is approximately 50% higher than for a 
singular lens because the average magnification is slightly higher. For lower redshift sources 
the limits on the core radius will be somewhat smaller, and for higher redshift sources the 
limits will be somewhat higher. 

Although a large core radius can be used to escape the lensing constraints, such an 
assumption leads to some contradictions. If the core radius is tuned to reduce the number 
of lenses by an order of magnitude, then we expect all observed lenses to be marginal lenses. 
Marginal lenses generally have visible central or "odd images," and their absence can be 
used to set limits on the core radius (Wallington & Narayan 1993, Kassiola & Kovner 1993). 
In the observed sample of large separation quasar lenses we do not see any central images, 
suggesting that the core radii of clusters must be near singular. An alternate, and logically 
dangerous, argument is to reverse the direction of argument and note that when we normalize 
the cosmogonic model based on other observations (COBF, cluster number counts, and 
correlation functions), we only produce the observed numbers of lenses if the clusters are 
assumed to be nearly singular. 

7 Conclusions 

The key to making a quantitative comparison between observations of large separation 
gravitational lenses and cosmological predictions is a good model for the selection effects 
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Fig. 9. -The integral separation distributions (left panel) and the differential lens redshift prob- 
ability distribution for 0957+561 (right panel) for various models with h = 0.5. The solid line 
(COBE CDM) shows the distributions for COBE normalized CDM (o-g = 1.4) and the solid line 
with points (lens CDM) shows the distributions for CDM normalized to produce 3 lenses larger 
than 5'.'0 (ag = 0.5). In the left panel, the number of lenses produced by the COBE normalized 
CDM model is divided by a factor of 10 to fit it on the same scale as the other models. The figure 
also shows a tilted model (n = 0.6, ag = 0.6, dashed lines), a T model (F = 0.25, ag = 0.6, dashed 
lines with points), an open model (fio = 0.45, ag = 0.8, heavy solid lines), and a fiat model with 
a cosmological constant (fio = 0.35, ag = 0.9, heavy solid lines with points) that produce approxi- 
mately 3 lenses larger than 5'.'0 and are consistent with COBE. The tilted model and the F model 
give nearly identical results. The histogram in the left panel shows the integral separation distribu- 
tion of the known lens Q0957-|-561 and the two best candidates Q2345-|-007 and Q1635-|-667. The 
vertical line in the right panel marks the known redshift of the lens in 0957-|-561. 

that determine whether the lenses are found by quasar surveys. Based on the simple picture 
that quasar surveys must have optical magnitude limits, we built a selection function model 
that explains the magnitude distribution of wide separation pairs in the HB93 catalog. One 
important implication of the model is that quasar surveys are a very inefficient method of 
finding wide separation quasar pairs. 

We use the Press-Schechter (1974) approximation to compute the expected number of 
lenses in different cosmogonic scenarios, an approach pioneered by Narayan & White (1988). 
When we compare the PS results to the Wambsganss et al. (1994) numerical simulations we 
find that the overall agreement is good. For example, the difference in the expected number 
of lenses with separations larger than lO'.'O for a source at reshift Zs = 2 is only 20%. Such 
small differences have negligible effects on estimates of CTs, the rms fiuctuation over a 8h~^ 
Mpc top-hat window function, because the number of lenses depends exponentially on ag. 
Some differences between the PS models and the Wambsganss et al. (1994) simulations are 
clearly traceable to differences in the models. The lens redshift distributions differ because 
we use a singular isothermal sphere for the lens model, while the simulations have mass 
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distributions with finite core radii (~ 10h~^ kpc) set by the resolution of the calculation. 
When we compute the redshift distribution with a f0-20/i~^ kpc core radius for the lenses 
we reproduce the numerical simulations. For large separation lenses the PS results predict 
many more wide separation lenses, and the deviations between the two results begin with 
separations of order 40'.'0. This is approximately the image separation scale corresponding 
to perturbations comparable to the outer scale of the numerical simulations [5h~^ Mpc). 
Other differences between the models are not easily assigned to either approach. 

Systematic errors in the PS approach, either in estimating the critical overdensity for 
collapse, Sc{z)^ or the virial radius of the collapsed perturbation do not strongly modify the 
conclusions. The estimates of erg vary linearly with rescaling either 6c or cr/ro, so rescaling 
either of these relations by 50% rescales the estimate of erg by 50% at erg — 0.5. Changing 
the expected number of lenses by a factor of two either by altering the selection function 
model or the structure of the lenses changes the estimate of CTs by 0.1 for erg — 0.5. Adding 
ellipticity to the lens model leads to a small change in the expected number of lenses, and 
many (30% to 50%) of the lenses detected by quasar surveys will consist of two images of a 
four image lens. Since the actual number of large separation lenses is small and uncertain 
(at least one, possibly three) most of theses systematic errors change the estimated value of 
(Ts by amounts comparable to the intrinsic statistical uncertainties. 

The addition of a large core radius to the density distributions of clusters and groups can 
strongly affect the conclusions. A core radius larger than 15h~^cr^ kpc reduces the expected 
number of lenses by a factor of two, and a core radius larger than 19h~^crl reduces the 
expected number of lenses by a factor of ten. Simulations of dark matter halos (eg. Dubinski 
& Carlberg 1991), galaxy scale lenses (Wallington & Narayan 1993, Kassiola & Kovner 1993), 
and the models of giant arcs (see Soucail & Mellier 1994) all suggest that potentials tend to 
have compact core radii. Moreover there is a serious fine tuning problem in giving groups 
and clusters core radii large enough to allow consistency with COBE normalized CDM, while 
producing only a few lenses and making the observed lenses have undetectable odd images 
in the lens cores. Neither the PS simulations nor the numerical simulations are capable of 
addressing this issue in any detail, but a larger sample of lenses will strongly constrain the 
core radii. 

We determined the values of erg required to produce the observed number of wide sepa- 
ration lenses in standard CDM, tilted CDM, low Oq cosmologies, and fiat cosmologies with 
a cosmological constant. Figure 9 summarizes the distributions of lens separations and the 
probability distribution for the lens redshift for standard CDM and various models normal- 
ized to fit both COBF and the lens data. The general results agree with other estimates 
of (Ts in the various models. For standard CDM models we find erg ~ 0.45 it 0.2 {2cr) for 
h = 0.5 similar to estimates from correlation functions (Maddox et al. 1990, Picard 1991, 
Vogeley et al. 1992, Loveday et al. 1992) and cluster abundances (Peebles et al. 1989, Frenk 
et al. 1990, Bahcall & Cen 1992, 1993), and in strong disagreement with the normalization 
erg = (2.8 ± 0.2)h needed to fit the COBF observations (Smoot et al. 1992, Wright et al. 
1994, Gorski et al. 1994). As previously noted by Narayan & White (1988), Cen et al. 
(1994), and Wambsganss et al. (1994), COBF normalized CDM overpredicts the number of 
lenses by more than an order of magnitude. It also predicts a much fiatter distribution of 
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image separations than is observed, and higher average lens redshifts (see Figure 9). Tilted 
models with a primordial spectrum oc k^ with 0.3 ^ n <^ 0.7 (h = 0.5) are consistent with 
the COBE estimates for erg at those values of the exponent n. The limits on the exponent n 
from the COBE data are n = 1.2 ± 0.3 (Wright et al. 1994, Gorski et al. 1994), so the over- 
all consistency is poor. Changes in the Hubble constant have negligible effects on the lens 
models, with the best fit value for erg decreasing by 0.1 when we increase h from 0.5 to 1.0. 
Treated purely as a fitting function in the power spectrum, effective Hubble constants (the 
r parameter in the power spectrum) of 0.15 ^ F ^ 0.30 are consistent with the COBE data 
and the observed number of lenses. All the variant models that produce the observed number 
of lenses and are consistent with COBE have steeper distributions of lens separations and 
lower average lens redshifts than standard CDM (see Figure 9). Since lensing determines erg 
of the mass distribution, the fact that it agrees with estimates of erg from luminous objects 
suggests that there is little or no bias on cluster scales. We did not examine cold + hot 
models, although they should show better agreement because they reduce the numbers of 
groups and clusters (eg. Nolthenius et al. 1994). 

Models with low matter densities whether fiat with a cosmological constant or open 
models are consistent with COBE and the lens observations if 0.3 ^ Oq ^ 0.5 for h = 0.5. 
Unlike galaxy scale lenses where evolution appears to be unimportant (Mao 1991, Mao & 
Kochanek 1991, Rix et al. 1994), evolution matters far more than the extra volume of 
the low matter density universes. The same number of lenses is produced by higher values 
of (Ts largely because the mass of a collapsed perturbation for a fixed comoving scale Tq 
decreases rapidly when Oq is reduced. The high cosmological constant models are strongly 
ruled out by gravitational lensing (Kochanek 1993b, Maoz & Rix 1993, Kochanek 1994) with 
a current two standard deviation limit of Aq ^ 0.6. The low value of the COBE quadrupole 
compared to the higher £ fiuctuations sets a weaker limit that Aq ^ 0.8 (Sugiyama & Silk 
1994). Moreover, if the cosmological constant is invoked to solve the age problem produced 
by a high value of the Hubble constant, the value of Aq is forced to be larger than 0.85 
exacerbating the confiict with the other limits. 

In the models normalized to produce the observed number of lenses, the average com- 
pleteness of the lens sample is 20%. In the best fitting models, 90% of the lenses are smaller 
than 20'.'0. The quasar surveys miss the lenses because of their limited average dynamic 
range, and the existing lens surveys are poorly designed for finding wide separation lenses. 
Fens surveys examine regions 5" to 10" in radius around each quasar for lensed images using 
one or two color photometry followed by spectroscopy. The survey region is limited because 
the background of galactic stars is already several times the number of lenses within these 
small regions. A wide separation lens survey is equivalent to a new, deep quasar survey 
centered on a known quasar, with the same selection problems and contamination problems 
of quasar surveys. If the goal is restricted to finding lensed images of a known quasar, the 
optimal solution for surveying large areas around the quasar is to use intermediate-width 
filters (200A wide) bracketing the strongest available emission line of the quasar. Since the 
equivalent width of a strong quasar emission line is a large fraction of the filter bandwidth, 
there should be very few objects other than quasars at the same redshift capable of matching 
the color of the quasar. 
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